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Abstract 

Feral honey bee populations have been reported to be in decline due to the spread of Varroa destructor, an ectoparasitic 
mite that when left uncontrolled leads to virus build-up and colony death. While pests and diseases are known causes of 
large-scale managed honey bee colony losses, no studies to date have considered the wider pathogen burden in feral 
colonies, primarily due to the difficulty in locating and sampling colonies, which often nest in inaccessible locations such as 
church spires and tree tops. In addition, little is known about the provenance of feral colonies and whether they represent a 
reservoir of Varroa tolerant material that could be used in apiculture. Samples of forager bees were collected from paired 
feral and managed honey bee colonies and screened for the presence of ten honey bee pathogens and pests using qPCR. 
Prevalence and quantity was similar between the two groups for the majority of pathogens, however feral honey bees 
contained a significantly higher level of deformed wing virus than managed honey bee colonies. An assessment of the 
honey bee race was completed for each colony using three measures of wing venation. There were no apparent differences 
in wing morphometry between feral and managed colonies, suggesting feral colonies could simply be escapees from the 
managed population. Interestingly, managed honey bee colonies not treated for Varroa showed similar, potentially lethal 
levels of deformed wing virus to that of feral colonies. The potential for such findings to explain the large fall in the feral 
population and the wider context of the importance of feral colonies as potential pathogen reservoirs is discussed. 
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Introduction 

The feral honey bee colonies of the UK were thought to have 
been severely reduced by the arrival of the Varroa mite in 1992 
[1,2], which also caused severe losses in the managed population 
[3]. However, anecdotal reports have spoken of a recent 
resurgence in the feral honey bee population, and have suggested 
this could be due to the evolution of an adapted host parasite 
relationship [3-5]. 

Varroa destructor is an ectoparasitic mite that has both direct 
and indirect impacts on honey bee health. The mite causes direct 
damage to the developing honey bee larvae and pupae by sucking 
their haemolymph and reducing their hatching weight [5]. Bees 
parasitized in this way usually begin foraging earlier and have a 
significantly reduced lifespan which may be due to decreased 
learning abilities, impaired ability to navigate and consequently a 
lower probability of returning to the colony [5] . Indirect effects of 
V. destructor are termed varroosis, whereby the Varroa mite acts 
as a vector for a variety of honey bee viruses, most notably 
deformed wing virus (DWV) [6] . 

Before the occurrence of Varroa mites, bee viruses were 
generally considered a minor problem to honey bee health [5], 
Recently, however, Genersch (20 10) found DWV and ABPV to be 
significantly related to German winter colony loss, while Highfield 
et al (2009) attributed 67% of overwintering colony loss in 



Devonshire to DWV [7,8]. Indeed varroosis is now considered to 
be the most destructive disease of honey bees worldwide [5,6] and 
a major cause of winter colony loss [9,10]. 

Beekeepers control Varroa levels in colonies using synthetic 
acarides, organic acids, essential oils and a wide variety of 
management techniques, which has helped to improve survival 
rates [5,11]. It was assumed, given the early large scale losses of 
untreated managed populations that feral, and thus untreated 
colonies, must quickly succumb to varroosis, although no evidence 
was ever presented to support this theory. However, it has been 
also suggested that sufficient time has passed since the first 
exposure to Varroa mite infestation to allow selection pressure to 
act on bee populations, and that feral honey bee populations are 
starting to rebound [4, 1 2] . Indeed, shorter-term selective breeding 
of managed colonies for 'Varroa resistance' has been shown to 
lower Varroa numbers in some colonies [4,13,14]. In the Arnot 
Forest, in the US, feral colonies were found to be at the same 
density in 2007 as in pre Varroa times in 1978 [13,15]. Untreated 
colonies in France and Sweden were also used for a 'live and let 
die experiment', where colonies that survive without Varroa 
treatment are subsequently selected for honey production in an 
attempt to create a race that is both Varroa tolerant and 
economically attractive [13,14,16]. There is also a strong body 
of evidence from laboratory and field studies to show rapid 
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evolutionary responses to selection pressures in other insect-pest 
interactions, such as in the case of Drosophila melanogaster 
populations and the parasitoid Asobara tabida [17,18]. 

If found to be coping with varroosis in the absence of active 
management, feral honey bee colonies could act as important 
genetic stocks from which to improve breeding efforts for mite 
tolerant managed honey bees [5,19]. Alternatively, feral nests 
could simply represent escaped swarms from managed colonies 
that could present a risk to the managed population by harbouring 
disease agents and re-infecting managed stocks [20,21]. 

Controlling communicable disease in managed honey bee 
populations is a challenge, given honey bees can move disease 
agents over great distances. Adult bees can be used to infer the 
infection state of a colony, allowing the disease state of a colony to 
be determined without the need for a destructive sample of brood 
[22]. This is the first time these methods have been used to 
measure the pathogen burden of the UK's feral honey bee 
population, and the first study to compare feral pathogen levels 
with those of local managed colonies. 

Methodology 

Site selection 

Feral honey bee colonies were located by engaging the 
beekeeping community and the general public using several 
methods; (1) emails to beekeeping associations (both the main 
BBKA secretary, but also secretaries of regional beekeeping 
associations); (2) notes on applicable internet forums such the 
natural beekeeping forum; and (3) an article in BeeCraft, a popular 
monthly beekeeping magazine [23]. Respondents were able to 
report their colony by email, letter, or using a bespoke website 
(www.honeybeeproject.co.uk). 

Locations of feral colonies were selected based on a good history 
of activity at the nest site (1 year minimum) thus avoiding the 
inclusion of new swarms with no history of survival. Sites that were 
inaccessible were not selected for safety reasons. The managed 
apiary nearest the feral site was identified using a national 
beekeeping register called BeeBase (see www.nationalbeeunit.com) 
and samples of adult honey bees were collected from feral and 
managed colonies on the same day in Spring 2010. At the same 
time each beekeeper was asked whether they actively controlled 
Varroa. Unfortunately it was impossible to carry out Varroa 
screening at the colony site as most cavities were inaccessible and 
Varroa mites were only seldom seen on adult workers collected at 
the colony entrance. 

All sites surveyed were on private property. Permission was 
sought from the individual land owners for house and garden sites, 
and the National Trust for parks and estates. Field studies did not 
involve protected or endangered species. 

Nucleic acid extraction 

Approximately 60 foraging A. mellifera adults were collected 
from each colony and stored for use in 100% ethanol at — 70°C 
[24] . Twenty-four bees from each of the 34 paired colonies were 
selected. Whole bees were washed in molecular grade water, and 
individually disrupted with 2.3 mm silica beads in bead beater at 
8000 g for 30 seconds (Precellys). Total DNA and RNA was 
extracted from each worker bee using a 10% Chelex solution with 
TE buffer. After disruption, 800 |xl of 10% Chelex solution was 
added to each crushed bee residue. The solution was heated to 
95°C for 5 minutes then centrifuged at 8000 g for a further 5 
minutes. The upper aqueous layer was removed (200 ul) before 
centrifugation at 8000 g for 5 minutes, with the removal of 150 u.1 
of the upper aqueous DNA. Finally, 20 ul of extract from each 



individual bee was pooled per colony [8]. Colony extractions were 
divided and purified to optimise RNA [25] and DNA recovery 
[26]. Briefly, 300 ul of extract was mixed with 300 ul of 24:1 
chloroform: IAA solution and spun at 8000 g for 10 minutes. For 
RNA, 100 ul of the upper aqueous layer was transferred into a 
fresh tube containing 1 00 ul of 4 M LiCl and samples were mixed 
well and left overnight. For DNA 100 u.1 of the upper aqueous 
layer was transferred into a fresh tube containing 50 |J.l of 5 M 
NaCl and 100 u.1 isopropanol. Next both DNA and RNA samples 
were vortexed and centrifuged for 10 minutes 8000 g. The 
remaining salt and ethanol was decanted and the pellet was 
washed with 500 u.1 of 70% ethanol by spinning for 4 minutes. The 
ethanol was decanted and the pellet dried in a heated vacuum for 
5 minutes at medium heat. Dried pellets were re-suspended in 
150 ul of lx TE buffer and frozen at -20°C until required. 
Extraction blanks were created during this process using the same 
reagents. 

Parasite testing 

Colony extracts were tested for the presence of Acarapis mites, 
black queen cell virus (BQCV), chronic bee paralysis virus 
(CBPV), deformed wing virus (DWV), Nosema apis, Nosema 
ceranae, Melissococcus plutonius, Paenibacillus larvae, Slow 
paralysis virus (SPV) using established protocols (Table 1). Real- 
time reactions were set-up using TaqMan chemistry using PCR 
core-reagent kits (Applied Biosystems, Branchburg, New Jersey, 
USA), according to the manufacturer's protocols. Reactions (25 u.1) 
were comprised of 2.5 u.1 of buffer (Buffer A), 5.5 ul MgCl 2 
(25 nM), 2 ul dNTP, 1 u.1 of forward and reverse primers, 0.5 |xl of 
probe, 0.125 ul Taq polymerase, 0.125 ul MMLV (RT reactions 
only) and 5 ul of DNA extract (parasites) or RNA (viruses). All 
Taqman probes were specific and covalently labelled with a 
reported dye (FAM) at the 5' end and with a quencher dye 
(TAMRA) at the 3' end (Table 1). Samples were run in triplicate 
reactions with template positive and extraction blanks. 

Reactions were run on an ABI Prism 7900HT (Applied 
Biosystems) with real-time data collection. Reverse transcription 
was performed at 48°C for 30 minutes, followed by denaturing 
and enzyme activation at 95°C for 10 minutes. This was followed 
by 40 cycles of denaturing at 95°C for 15 seconds and a combined 
annealing and extension step for 60 seconds at 60°C. Fluorescence 
values, amplification plots and threshold cycle (C t ) values were 
calculated using SDS 2.2 (Applied Biosystems). 

Quantification of PCR results 

The relative amount of DWV, BQCV, N. apis and N. ceranae 
were analysed using the comparative C, method [27] using the 
parasite specific assay as the target and Elongation Factor 1 (EF1) 
as the reference assay (Table 1). The use of this gene as an internal 
control is established [28,29]. A survey of control genes, suggested 
that variation among the samples using four different references 
genes for normalization, including EF1, was rather small and did 
not drastically change the target gene expression profiles between 
samples [30] . PCR efficiencies between target and reference assays 
was compared by diluting a known positive sample 1:10 through 6 
levels. Efficiencies above 90% with an R2 above 0.98 indicated 
data [31]. 

Wing morphology 

Wing morphometry data was gathered for all colonies where 
more than 50 useable right forewings were available. Wings were 
removed from the bee and placed under a glass on an Epson 
Perfection V300 Photo scanner. Images were scanned at 4800dpi 
resolution using positive film strip mode. DrawWing software 
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Table 1. Primers used in this study. 



Target 


Primer name 


Sequence (5 -3 ) 


Acarapis spp. 5 


Acarapis F1 


GCCATAAGACATCACTATTCT 




Acarapis Rl 


TCATTTAAACTTCATGATACTCTCAATCAG 




Acarapis 1 


TGCGCAATGCAACTAGTCCTCTAAAGACTAGTTTC 


Black queen cell virus 1 


BQCV8I95F 


GGTGCGGGAGATGATATGGA 




BQCV 8265R 


GCCGTCTGAGATGCATGAATAC 




BQCV8217T 


TTTCCATCTTTATCGGTACGCCGCC 


Chronic bee paralysis virus 1 


CBPV 304F 


TCTGGCTCTGTCTTCGCAAA 




CBPV 371R 


GATACCGTCGTCACCCTCATG 




CBPV 325T 


TGCCCACCAATAGTTGGCAGTCTGC 


Deformed wing virus 1 


DWV 958F 


CCTGGACAAGGTCTCGGTAGAA 




DWV9711R 


ATTCAGGACCCCACCCAAAT 




DWV 96271 


CATGCTCGAGGATTGGGTCGTCGT 


Elongation factor I 2 (Internal control) 


FF1 F 


CTGGTACCTCTCAGGCTGATTGT 




FF1 R 


GCATGCTCACGAGTTTGTCCATTCT 




FF1 T (TAMRA) 


TGCTTCGAACTCTCTCCAGTACCAGCAGCAACA 


Nosema apis 5 


N apis Fl 


ATTTACACACCAGGTTGATTCTGC 




N apis Rl 


TGAGCAGTCCATCTTTCAGTACATAGT 




N apis MGB 


TGACGTAGACGCTATTC 


Nosema ceranae 5 


Nosema cl 83G F 


TTG AGA GAA CGG TTT TTT GTT TGA G 




Nosema cl 974 R 


TTC CTA CAC TGA TTG TGT CTG TCT TTA A 




Nosema cl 865 T 


ATA ATA GTG GTG CAT GGC CGT TTT CAA TGG 


Melissococcus plutonius 3 


EFBF 


TGT TGT TAG AGA AGA ATA GGG GAA 




EFB Rev2 


CGT GGC TTT CTG GTT AGA 




EFBProbe 


AGA GTA ACT GTT TTC CTC GTG ACG GT 


Paenibacillus larvae 5 


PI_R24_468F 


TCCCCGAGCCTTACCTTTGT 




PI_R24_538R 


ACCTACGAACTTGACGCTGTCCT 




PI_R24_489T 


TGCTCATACCCGGTCAGGGATTCGA 


Slow paralysis virus 4 


SPV 8383F 


TGATTGGACTCGGCTTGCTA 




SPV 8456R 


CAAAATTTGCATAATCCCCAGTT 




SPV 84071 (TAMRA) 


CCTGCATGAGGTGGGAGACAACATTG 


Sacbrood virus 7 


SBV3IIF 


AAG TTG GAG GCG CGy AAT TG 




SBV 380R 


CAA ATG TCT TCT TAC dAG AGG yAA GGA TTG 




SBV33I1 


CGG AGT GGA AAG AT 


The 5'-terminal reporter dye for each TaqMan 


probe was 6-carboxyfluorescin 


TAM) and the 3' quencher was tetra-methylcarboxyrhodamine (TAMRA) or Minor groove 



binding (MGB) as indicated. 
'-[•«]. 

2 - m. 

3 - [48]. 
"-[50]. 
5 -[51]. 



doi:10.1371/journal.pone.0105164.t001 

Table 2. Results from testing nucleic acid preparations for parasites with fewer than five positive observations. 





Disease 


Feral Positives 


Managed Positives 


Acarapis spp. 


1 


4 


CBPV 


3 


4 


M. plutonius 


0 


1 


P. larvae 


0 


0 


SBV 


1 


0 


SPV 


0 


0 



doi:1 0.1 371 /journal.pone.01 051 64.t002 
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Figure 1. The Restricted Maximum Likelihood model estimates for the four most commonly found pathogens. Predictions are on the 
log scale with 95% confidence intervals. * denotes a significant different between paired managed (m) and feral (f) colonies. Analyses are done 
separately for each pathogen and not between pathogens. 
doi:10.1371/journal.pone.0105164.g001 



version 0.45 was used as the best example of modern wing 
morphometry, to record the cubital, hantel and discoidal shift 
index [32,33]. Landmarks were placed by eye were DrawWing 
failed to correctly identify venation junctions on wings with slight 
damage to the tips. DrawWing outputs were transferred to 
Morphplot version 2.2 [34] and the average cubital, hantel and 
discoidal shift index calculated for each colony. 

Statistical analysis 

Data for the four most commonly found parasites were analysed 
by Restricted Maximum Likelihood (REML) to account for the 
paired structure of the data using GenStat 14.1 [35]. The pairs 
were included in the model as a random effect whilst the treatment 
of interest (managed vs feral colonies) was included as a fixed 
effect. Further, the data were log-transformed to correct for right 
skew. 

Log DWV levels were compared between three groups: feral 
(untreated), managed (treated) and managed (untreated) by 
analysis of variance. Individual treatments were compared using 
post-hoc test and a Tukey adjustment for multiple comparisons. 
This analysis was carried out in R version 3.0.0 [36]. 

Cubital index, discoidal shift angle and hantel index was 
analysed for feral and managed colonies. These index values were 
log-transformed to account for a right skew, and analysed by 
multivariate analysis of variance, using R version 3.0.0. 



Results 

Paired samples 

In total, 1 00 reports of feral colonies were received and 60 were 
visited in Spring 2010. Of those visited, 34 feral colonies were 
sufficiently accessible to collect foraging bees and each was 
successfully paired with a nearby site containing a managed honey 
bee colony. Pairs were an average of 1 .4 km apart, with a 
maximum of 10.9 km. 

Parasite testing 

C t values for positive controls and buffer blanks were checked 
initially, to ensure correct functioning of the master mix and an 
absence of contamination. All colony samples tested negative for 
SPV and P. larvae, M. plutonius, SBV, Acarapis spp. and CBPV 
were detected at low prevalence (Table 2). 

All colonies were positive for DWV, BQCV, N. apis and N. 
ceranae. Similar PCR efficiencies (between 92% and 107%) were 
achieved for target assays and 96% for the reference assay allowing 
quantitative interpretation of the commonly detected pathogens 
[37]. There was no significant difference in the amount of N. apis 
(F= 1.70, d.f= 1, 33, p = 0.20,), N. ceranae (F = 0.52, d.f. = 1,33, 
p = 0.48) or BQCV (F = 1.11, d.f. = 1,33, p = 0.30) between feral 
and managed colonies. However, the amount of DWV was 
significandy different between managed and feral colonies 
(F = 6.41, d.f= 1,33, p = 0.016) (Fig. 1, 2 and 3). 
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Figure 2. The proportion of DWV between pairs of either feral/managed or feral/untreated managed colony pairs. The pie charts are 
split into two groups: the red fill indicates feral colony DWV levels alongside the blue fill for managed colonies where no Varroa treatment was used. 
The black fill indicates feral colony DWV levels alongside the white fill of Varroa treated managed colonies. 
doi:1 0.1 371 /journal.pone.01 051 64.g002 



In total, usable right wings were obtained in sufficient number 
from 56 colonies (28 feral and 28 managed) to obtain wing 
morphometry metrics. There was no significant difference 
between wing morphometric indices for managed (n = 28) and 
feral colonies (n = 28) (approx. F = 0.43; df = 6, 104; p = 0.86, 
Fig. 4). 

Discussion 

We present novel data to describe the parasite burden of feral 
honey bee colonies and report that three of the major parasites 
show similar intensity in both managed and feral adult honey bee 
populations. Crucially, the levels of DWV were far higher (Fig. 1) 
in feral colonies compared to managed colonies, which could 
reflect the absence of Varroa control in feral colonies. 

Interestingly managed colonies not treated for Varroa contained 
similar virus levels to feral colonies (Fig. 3). DWV is the most 
prevalent known virus in managed honey bee colonies in Europe 
and an increasing proportion of viruliferous Varroa mites has been 
linked to reduced colony survival [10]. Colonies with high levels of 
DWV show evidence of a scattered brood nest, crippled bees, loss 
of coordinated social behaviour such as hygienic behaviour, queen 
attendance and a rapid decline in the colony's bee population 
[5,6,10,38]. The extremely high values of DWV found in feral and 
untreated managed colonies would be expected to lead to colony 
mortality [8], however it is not clear whether feral colonies have 
novel mechanisms to resist such high levels of DWV. Whilst it is 
not possible to rule out increased Varroa tolerance or a balanced 
host parasite relationship in feral nests [3-5], the development of 
alternative coping strategies to mitigate the high DWV levels 
detected seems unlikely. 



For the first time, we also present results that suggest the 
venation of the forewings of feral honey bees is not distinguishable 
from managed honey bees using three different morphometric 
measures (Fig. 4). Both the native honey bee of the British Isles 
(Apis mellifera mellifera), and the dominant European races (A, m. 
carnica, A. m. iberica and A. m. ligustica;[39]) all have well 
characterised morphometry. Whilst wing morphometry is often 
used to determine honey bee race, and has been criticised for 
being over reliant on the wing section of the genome [40] , it gives 
us a useful first assessment of the apparent genetic similarity 
between feral and managed colonies [41]. Such morphological 
similarities suggest that the feral population may simply be a 
consequence of escapees from the managed population. Anecdotal 
increases in feral populations could simply reflect a combination of 
a recent increase in the popularity of beekeeping leading to a 
higher number of novice beekeepers who are likely to allow their 
colonies to swarm, and an increasing proportion of let-alone 
beekeepers who actively encourage natural honey bee behaviours 
like swarming [42] . A workable level of Varroa tolerance is keenly 
sought by UK beekeepers, but abandoning Varroa treatment 
without ensuring colonies have evolved a natural resistance to the 
mite, or without virus-free Varroa, could leave colonies danger- 
ously exposed, particularly in areas of high beekeeping density 
[7,43]. Indeed the geographic proximity, and lack of large, remote 
and untreated honey bee populations may prohibit meaningful 
breeding programs for Varroa resistance in England and Wales. 
Instead, a more gradual approach of selective breeding for Varroa 
tolerance is likely to lead to the improvements in resilience that 
global apiculture requires [44]. 

More BQCV was found in feral colonies compared to managed 
colonies, but this difference did not reach statistical significance 
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Figure 3. The effect of Varroa treatment on managed and feral colony log DWV levels separated by treatment. Blue indicates feral 
colonies untreated for Varroa. Red indicates managed colonies undergoing standard Varroa treatment (i.e. dosing with Varroacide one to two times 
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(Fig. 1). V. destructor has been reported to vector this virus [37], so 
it is not inconceivable that this slight increase is also linked to mite 
parasitism. Nosemosis C is an emerging disease in Europe caused 
by N. ceranae [45], yet is found to be well established in feral 
colonies. This result could suggest a degree of parasite perturba- 
tion between feral and managed populations [46]. One can 
hypothesise that feral honey bees are exposed to fewer stressors in 
the form of beekeeper manipulation e.g. direct damage to comb 
and propolis, death of bees during beekeeper activity, cross 
contamination between hives, honey removal, pollen harvesting 
etc. [16]. In this study only one colony was positive for the 
causative agent of European foulbrood, a disease sometimes linked 
to stress [47], and this was a managed colony. The prevalence of 
EFB is low in the UK [48], so the small sample size makes it 
impossible to draw any meaningful conclusions about the true 
level of foulbrood infection in the UK's feral colonies. However, 
we can conclude that the sampled feral population contained 
many of the parasites commonly found in managed populations. 
Therefore, our study provides no evidence that feral-nests reduce 
parasite load compared to managed nests. However, our results do 
not address the consequences of the measured pathogen burden in 
feral nests, and it remains possible that feral colonies are more able 
to cope with the observed pathogen load than their managed 
counterparts. 

Given the novel observations that (i) feral colonies contain 
crippling high levels of DWV; (ii) managed and feral populations 
appear similar using three different measures of wing morphom- 
etry and (iii) feral and pathogen populations share even recently 



emerged parasites, it seems likely that the invasion of the Varroa 
mite and the increase in prevalence of its concomitant viruses may 
indeed explain the loss of feral honey bee colonies. Despite 
showing high levels of DWV in feral colonies, we cannot 
categorically link this to an increase in feral colony mortality. 
Future studies could concentrate on understanding whether our 
observations of high DWV titre result in colony mortality or 
whether feral populations have behavioural adaptations, such as 
increased swarming, to tolerate levels of DWV that would be 
detrimental to a managed colony. Finally, future work could use 
microsatellite markers to categorically explore the relatedness of 
feral and managed honey bee populations. 
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